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ABSTRACT 

Here  we  solve  the  equation  of  Magneto  hydrodynamic  boundary  layer  flows  past  a moving  continuous  flat  surface 
in  the  presence  of  transverse  magnetic  field.  Here  we  check  effect  of  magnetic  field  on  velocity  of  MHD  fluid.  Solve  the 
governing  nonlinear  equations  with  use  their  associated  boundary  conditions  with  blue  method.  The  beauty  of  this  method 
is  we  can  solve  nonlinear  problem  directly,  without  convert  in  linear  form.  Numerical  solution  and  graphical  presentation 
of  flow  problem  is  also  given. 

NOMENCLATURE 

MHD  - Magneto  hydrodynamics 

U,  V - Velocity  components  in  X , Y directions  respectively 

U,W  - The  mainstream  velocity  components 

p - Field  density 

\j/  - Stream  function 

<7  - Electrical  conductivity 

By  (a)  -magnetic  field  strength 

Vm  - Magnetic  viscosity 

S - Magnetic  field  strength 
T]  - The  single  independent  variable 

KEYWORDS:  Boundary  Layers,  Blue  Method,  Linear  Equations,  Magnetic  Field,  Third  Order  Nonlinear  Differential 
Equation 

1.  INTRODUCTION 

Due  to  wide  application  in  several  technical  and  industrial  processes,  the  boundary  layer  flow  over  continuous 
moving  surface  have  attracted  many  researcher  in  many  branches  of  engineering.  MHD  is  the  motion  of  an  electrically 
conducting  fluid  in  the  presence  of  a magnetic  field.  Due  to  the  motion  of  an  electrically  conducting  fluid  in  a magnetic 
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field  the  electrical  currents  produces  their  own  magnetic  field  and  these  modify  the  original  magnetic  field.  This  induced 
current  interacts  with  the  magnetic  field  to  produce  electromagnetic  forces  effects  the  original  motion.  Thus  the  two 
important  basic  effects  of  Magneto  hydrodynamics  are  (i)  the  motion  of  the  fluid  affects  the  magnetic  field  and  (ii)  the 
magnetic  field  affects  the  motion  of  the  fluid.  MHD  flow  over  moving  surfaces  emerges  in  a large  variety  of  industrial  and 
technological  applications.  It  has  been  investigated  by  many  researchers.  [1]  Wu  (1973)  has  studied  the  effects  of  suction 
or  injection  on  MHD  boundary  layer  flow.  [2]  Takhar  et.al.  (1987)  studied  a MHD  asymmetric  flow  over  a semi-infinite 
moving  surface  [3]  Mahapatra  and  Gupta  (2001)  studied  a steady  two-  dimensional  stagnation-point  flow  by  a uniform 
transverse  magnetic  field.  [4]  Jean-David  Hoernel  (2008)  has  been  investigated  the  steady  laminar  incompressible 
boundary  layer  governing  MHD  flow.  [5]  Rajput  (2013)  studied  of  MHD  boundary  layer  flow.  [6]  Sharma  Vineetkumar 
(2011),  Numerical  solution  of  two  dimensional  MHD  stagnation-point  flow.  [7]  Patel  M.  (2013)  had  discussed  group 
theoretical  method.  [8]  Pandya  Jigisha  et  al.  (2010)  had  studied  foppl-hencky  equation  using  spline  method. 

The  governing  equation  is  nonlinear  differential  equations,  which  is  solved  by  using  spline  collocation  blue 
method.  In  this  way,  the  paper  has  been  organized  as  follows.  In  section  2,  basic  equations  for  given  problem  and  in 
Section  3,  spline  collocation  blue  method  and  approximate  solution  for  the  governing  equations.  The  results  and  discussion 
are  in  section  4. Conclusions  are  summarized  in  section  5. 

2.  BASIC  EQUATIONS 

[9]Srivastava  et  al  (1987),  the  basic  equations  governing  the  motion  of  two  dimensional,  steady  incompressible 
viscous  fluids  past  continuous  surface  in  the  presence  of  transverse  magnetic  field  can  be  written  as  follow: 

ay  c,^dy/ 

dy  dxdy  r)x  dy2  dy3  dy 


Where  S (x)  = 


(TB;(x) 

P 


With  the  boundary  conditions: 


dy/ 

dy 


■U(x) 


dy/ 

dx 


■■  -V(x)  at  y = 0 


(2.2) 


dy/ 

dy 


= 0 at  y — » 


OO 


(2.3) 


Where  X co-ordinate  measured  along  the  plate,  y co-ordinate  normal  to  the  plate,  U and  V are  the  velocity 
component  of  the  fluid  in  X and  Y directions.  V - magnetic  viscosity,  p - fluid  density,  <7  - electrical  conductivity, 
B (x)  -magnetic  field  strength  along  the  plat  wall  and  acting  perpendicular  to  it,  y/  - stream  function. 


To  study  the  physical  phenomena  of  flow  under  consideration,  solve  above  equations  (2.1),  (2.2),  (2.3).  Reduce 
the  partial  differential  equations  into  ordinary  differential  equations.  Patel  M.  [7]  reduces  the  partial  differential  equation 

into  ordinary  differential  equation.  For  that  they  employ  following  one  parameter  linear  group  transformation  F(  . 
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Here  F[  is  selected  as  (6) 


x = Aa'x,  y = Aaix 

Tl-.y  = Aa:ij7 , S (x)  = Aa>S(x) 

U = Aa5U,  V = Aa6V 


Where  OCi  '.i  = 1 6 and  A are  real  arbitrary  constant. 

Now  seek  relations  (X ' S such  that  the  basic  equations  (1)  along  with  the  boundary  conditions  (2)  and  (3)  will  be 
invariant  under  this  group  transformation  Fj . 

Invariant  conditions  demand  that  the  power  of  A in  each  term  of  transformed  equations  should  be  equal,  which 

yields 


2 a3  - 2 a2  - ax » a3  - 3 a2  = a4+a3-a2  (2.4) 

a3-a2=  a5  (2.5) 


a3-ai=  a()  (2.6) 

Solving  (2.4),  (2.5)  and  (2.6)  we  get 

df,  (X-,  cXi  . (Xc  , . (Xt 

= m(say),  — = 1 - m,  — = —2m,  — = 1 - 2m,  — = —m  (2.7) 

(X  (X  (X \ 


3.  QUARTIC  SPLINE  BLUE  METHOD 

For  three  points  boundary  value  problems  are 

Let  S{  (x)  be  quartic  spline  in  \xt  ] , X ] . Conditions  for  natural  splines  are 
i.  Sj  (x)  Almost  quartic  in  each  subinterval  [xj_l , X{  ] . 


ii.  Si(xi)  = yi  , for  i = 0,1,2, ,n. 

iii.  ,s’; (xj ) , s,  (,y),  s(.  (,y),  s (,y)  are  continuous  in  [x0 

iv.  s'(x0)  =s,m(xn)  = 0. 


Here  spline  third  derivative  must  be  linear  in  [x;_j , Xt . ] . 


So 


a"’(f)  = y-[(F  - x)  y.-i  + (*-  *>- 1 ) yi] 


(3.1) 
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Where  lll  — X!  — XM  and  S (x;  ) = y. 
Integrate  (1),  twice  with  respect  to  X . 


= i 


y a-.)3  : 

y'  - '-1  y'  - i 


+ cj  (xj  - x)  + di  (x  - xM). 


Where  use  (xM)  = yt_x  and  St  (x;)  = yi  in  (1). 


We  get  constants  Cj  and  d] 
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(x-xM). 


(3.2) 


Integrate  (3.2),  once  with  respect  to  x. 


*/(*)  = T- 
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(3.3) 


Take  A,.  (xM  ) = yM  we  get  constants  £'( 


Id  , h . 

Where  e,  = yM  - y yM  + - yM  • 


Substitute  in  (3.3),  we  get 


5;(x)  = 


h. 


(x(. -x)4  ...  (x  x_j  )4  , 

■y,--t  + J y,- 


24 


24 
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(3.4) 
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Here  st  (x,  ) = (jcf+)  and  for  equal  intervals  we  have, 

7 2 

yM ' ~ 2 y'i  + yt~ i = — (yi+i  + 4 3'/'  + .Vi”) 

6 

And  for  si(xi~)  = si+1(xf+)  and  for  equal  intervals  we  have, 


(3.5) 


yi-i-y,  =~{yi  + yi-i) 


h „ 

+-(3^  ) 


(3.6) 


3.1.  Solution  by  Using  Collocation  Method 

[8]  Srivastava  et  al  (1987),  the  basic  equations  governing  the  motion  of  two  dimensional,  steady  incompressible 
viscous  fluids  past  continuous  surface  in  the  presence  of  transverse  magnetic  field  can  be  written  in  non-linear  equation  as 
follow: 

f~Pf  +Jf"=0,  (3.1.1) 

Subject  to  boundary  conditions  as, 

/ (0)  = l,/(0)  = 1,/  (0.5)  = 0.  (3.1.2) 

To  obtain  the  spline  solution,  we  begin  with  a assume  function  f (j])  = —JJ~  + 7] + \which  satisfy  given 
boundary  conditions  (3.1.2).  To  find  the  solution  of  equation  (3.1.1)  along  with  boundary  conditions  (3.1.2).  First  we  use 
f (if)  = — : 77  -+77  + I and  (3.1.1)  in  (3.5)  and  /l  = 0.1  , we  gate  different  values  of  yt  for  ( = 1,2, 3, 4. 


To  find  the  final  solution  we  use  (3.6)  for  different  values  of  i = 1, 2, 3, 4 respectively 
We  get  the  equations  as  follows 

h , , h3  

y0-yi  = --\.yi+y0}+—\-yi  +3  y0] 

h , , h3  

yi-y2  = --ly2+y^+—i-y2  +3Ti  ] 

^ = -|[y3’  + y2  +3y2"] 

h , h3  

^3-3^4  =--[^4  +^3] + ^[-^4  +3  y3] 


(3.1.3) 


To  substitute  yf  and  yf  for/7  = 0.1  in  (3.1.3).  We  get  four  unknown  and  four  equations.  Solve  that  equations 
using  Matlab. 
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Table  1:  Solution  of  Problem  using  Numerical  Method 


Numerical  Solution  with  Spline 

n 

f(rj)  for  j3  = 0 

/ (7)  for  jB  = 1 

f(T])  for  J3  = 2 

0.0 

1.0000000000 

1.0000000000 

1.0000000000 

0.1 

1.0875475000 

1.0868558333 

1.0861641666 

0.2 

1.1515783333 

1.1494116666 

1.1472450000 

0.3 

1.1943391666 

1.1906141666 

1.1868891666 

0.4 

1.2181966666 

1.2133299999 

1.2084633333 

0.5 

1.2255974999 

1.2203058333 

1.2150141667 

Graphical  Solution  of  Given  Problem 


Figure  1:  Velocity  Profile  Versus  lj  for  Different  Values  of  fB 


Figure  2:  / (Tj)  Versus  7]  for  Different  Values  of  jB 

4.  RESULTS  AND  DISCUSSIONS 

From  Figure-(l),  it  shows  that,  there  is  a significant  impact  of  magnetic  field  on  the  velocity  profile  of  the  flow. 
Here  for  different  values  of  fB  , we  find  velocity  profile.  From  Figure  it  is  clear  that  Value  of  jB  increase  the  velocity  of 
fluid  decrease.  It  means  that  for  electrically  conducting  fluid  in  a magnetic  field,  magnetic  field  increase,  velocity  of  fluid 
flow  decrease.  Magnetic  field  is  transversely  proportional  to  velocity. 

5.  CONCLUSIONS 

Here  we  find  the  generalization  of  blue  method  for  third  order  problem  and  solved  the  problem  using  blue 
method.  The  beauty  of  this  method  is  no  need  to  convert  nonlinear  problem  into  linear,  we  can  solve  directly  in  nonlinear 
form.  This  shows  that  spline  method  also  gives  nearest  and  accurate  results.  This  shows  the  reliability  of  the  method.  Thus 
we  can  solve  such  type  of  problems  using  blue  method  without  convert  nonlinear  problem  into  linear  form. 
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